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Abstract: Recently, Kapustin, Willett and Yaakov have found, by using localization techniques, 
that vacuum expectation values of Wilson loops in ABJM theory can be calculated with a matrix 
model. We show that this matrix model is closely related to Chern-Simons theory on a lens 
space with a gauge supergroup. This theory has a topological string large dual, and this 
makes possible to solve the matrix model exactly in the large N expansion. In particular, we find 
the exact expression for the vacuum expectation value of a 1/6 BPS Wilson loop in the ABJM 
theory, as a function of the 't Hooft parameters, and in the planar limit. This expression gives an 
exact interpolating function between the weak and the strong coupling regimes. The behavior at 
strong coupling is in precise agreement with the prediction of the AdS string dual. We also give 
explicit results for the 1/2 BPS Wilson loop recently constructed by Drukker and Trancanelli. 
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1. Introduction 

Localization techniques in supersymmetric gauge theories have produced in recent years explicit 
expressions for a variety of correlation functions. In [1], they were used to prove the longstanding 
conjecture [2, 3] that the vcv of a 1/2 BPS Wilson loop in = 4 super Yang-Mills theory can 
be computed by a correlator in a Gaussian matrix integral. This gives the celebrated formula 



for the planar limit of the Wilson loop in the fundamental representation. In this formula, Ii is 
a modified Bessel function and A = g^N is the 't Hooft coupling. (1.1) gives an exact expression 
in A that interpolates between strong and weak coupling (see [4] for a review). 

More recently, Kapustin, Willett and Yaakov [5] have applied localization techniques to 
ABJM theories [6, 7]. These are superconformal field theories in three dimensions based on a 
U{Ni) X U{N2) Chern-Simons theory coupled to matter, and they have large N AdS4 duals. 
In [5], it was shown that the calculation of vevs of 1/6 BPS Wilson loops in these theories 
can be reduced to a calculation in a matrix model, and they verified that the Wilson loop vev 
reproduces the calculations of [8, 9, 10]. However, the matrix model they obtained is highly 
nontrivial, and the comparison with the gauge theory calculation was done by performing a 
perturbative expansion up to two loops. Some partial results on the planar limit of the matrix 
model have been obtained in [11]. 

In this note we point out that the matrix model of [5] can be solved by relating it to the 
Chern-Simons matrix models introduced in [12], and in particular to the lens space matrix model 
studied in detail in [13, 14]. In [13] it was shown that the lens space Chern-Simons matrix model 
is the large N dual of topological string theory on a certain class of local Calabi-Yau geometries, 
providing in this way a nontrivial generalization of the Gopakumar-Vafa duality [15]. This means 
in particular that the large N limit of the model can be studied by using standard techniques 
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in mirror symmetry. It turns out that the matrix model of [5] can be regarded as a supergroup 
version of the lens space matrix model studied in [12, 13, 14]. Since large N duals describe matrix 
models as well as their supergroup extensions, we can use topological string theory on a local 
Calabi-Yau geometry to obtain exact results in ABJM theory. 

In this paper we use the solution of [13, 14] to obtain an exact expression for the planar limit 
of the vev of the 1/6 BPS Wilson loop constructed in [8, 10, 9], as a function of the two 't Hooft 
parameters of the theory. This analytic expression is relatively complicated but it can be written 
down explicitly in terms of elliptic functions, sec for example (4.40), (4.43) for the case in which 
the two gauge groups have the same rank. It interpolates smoothly between weak and strong 
coupling, and at strong coupling it agrees precisely with the AdS theory prediction obtained in 



The 1/6 BPS Wilson loop operator that wc study in the planar limit involves only one of 
the gauge connections in the U{Ni) x U{N2) quiver, and it is not the most natural operator 
from the point of view of topological string theory. It has been recently shown in [16] that it 
is possible to construct a 1/2 BPS Wilson loop operator in the ABJM theory which localizes 
precisely to the natural Wilson loop operator for Chcrn-Simons theory on L(2, 1) (extended 
to a supergroup). The vev of such an operator can be calculated to all orders in 1/iV, as an 
exact function of the 't Hooft parameters, by combining the solution of [13, 14] with the matrix 
model-inspired techniques of [17, 18]. In the topological string dual, the vevs of these 1/2 BPS 
operators correspond to open topological string amplitudes. 

Our exact formulae for the vevs of Wilson loops in ABJM theory have the flavour of mirror 
symmetry. There are two sets of coordinates for the parameter space: the "bare" coordinates 
and the flat coordinates. The flat coordinates are identified with the 't Hooft parameters. They 
are computed by period integrals and they can be related to the bare coordinates through the 
mirror map. The vev of the Wilson loop is naturally expressed in terms of bare coordinates, and 
one has to invert the mirror map in order to re-express it in terms of 't Hooft parameters. In 
retrospect, one could say that the result of [3, 2, 1] for the 1/4 BPS Wilson loop in A/" = 4 super 
Yang-Mills theory is comparatively easier since in their case the relevant matrix model is the 
Gaussian one, with a simpler moduli space, and where the 't Hooft parameter is simply equal to 
the bare coordinate. 

The organization of this paper is as follows. In section 2 we review very briefly the matrix 
model obtained in [5] . In section 3 we present the solution of the lens space matrix model building 
on [13, 14]. In section 4 wc explain the relation between the two matrix models, we derive the 
exact results for the 1/6 BPS Wilson loop in the planar limit, and we study its behavior both at 
strong and at weak coupling. Finally, some conclusions are presented in section 5. 

2. Matrix model for the ABJM Wilson loop 

The ABJM theory is a quiver Chern-Simons-matter theory in three dimensions with gauge 
group U{N)k X U(N)_k and J\f = 6 supersymmetry. The Chern-Simons actions have couplings 
k and —k, respectively, and the theory contains four bosonic fields Cj, I = I,- - ,4, in the 
bifundamental representation of the gauge group. One can consider an extension [7] with a 
more general gauge group U{Ni)k x U{N2)-k- A family of Wilson loops in this theory has been 
constructed in [8, 9, 10], with the structure 



[8]. 
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where is the gauge connection in the U{Ni)k gauge group, df{(Ni) is the dimension of the 
represenation R of U{Ni), x^{s) is the parametrization of the loop, and M/ is a matrix deter- 
mined by supersymmetry. It can be chosen so that, if the geometry of the loop is a line or a 
circle, four real supercharges arc preserved. Therefore, wc will call (2.1) the 1/6 BPS Wilson 
loop. A similar construction exists for a loop based on the other gauge group, 



Wr = J^^^T^^R^ exp J l^iAf.x'' + '^\x\M'jCiCJ ) ds (2.2) 



where is the U{N2)-k gauge connection. The planar limit of the vev of (2.1) was computed in 
[8, 9, 10], for Ni = N2 = N, in the fundamental representation R= n, and in the weak coupling 
regime A <C 1, where 

N 



A = - (2.3) 



is the 't Hooft parameter. The result is 



{Wu) = l + -^X^ + 0{\'). (2.4) 

On the other hand, in the strong coupling regime A S> 1, the Wilson loop vev can be calculated 
by using the large N string dual, i.e. type IIA theory on AdS4 x [8, 9, 10]. This gives the 
prediction^ 

(Wu) ~ e-^. (2.5) 

As in the case of the 1/2 BPS Wilson loop in = 4 Yang-Mills theory, the exact answer for 
the planar limit of this vev should interpolate between the weak coupling behavior (2.4) and the 
strong coupling prediction of the large N string dual, (2.5). 

A crucial step in finding such an exact answer was taken in the paper [5]. It was shown 
there, through a beautiful application of the localization techniques used in [1], that the vev of 
(2.1) can be computed as a correlation function in a matrix model. This matrix model is defined 
by the partition function 

J L"; f3 cosh2(Mi^) 



i=l j=l 



where the coupling gg is related to the Chern-Simons coupling k of the ABJM theory as 



(2.6) 



= f. (2.7) 



One of the main results of [5] is that 



(Wr) = {Ttr (e'^O)ABJM , (2-8) 



^To be precise, the dual calculation is made by considering a fundamental string in AdS4, which is a 1/2 BPS 
object. In [8, 10] it was argued that the the strong coupling behavior obtained in this way should apply to the 
symmetric Wilson loop W^''™ defined in (2.12) below. It should also apply to the 1/2 BPS Wilson loop constructed 
in [16]. However, as we will show in this paper, the leading exponential behavior is common to all these Wilson 
loops. 
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i.e. the vev of the 1/6 BPS Wilson loop (2.1) can be obtained by calculating the vev of the 
matrix e^^ in the matrix model (2.6). This was explicitly checked in [5] by computing the vev 
in the r.h.s. of (2.8) in the matrix model, for the fundamental representation. Notice that the 
Wilson loop for the other gauge group, 

(Wr) = (Tr^ (e'^O)ABJM (2-9) 

is obtained from (2.8) simply by exchanging Ni ^ N2 and changing the sign of the coupling 
constant g —g. 

The Wilson loop (2.1) breaks the symmetry betwen the two gauge groups. Recently, a class 
of 1/2 BPS Wilson loops has been constructed in the ABJM theory [16] which treats the two 
gauge groups in a more symmetric way. These loops have a natural supergroup structure in 
which the quiver gauge group U (iVi) x U {N2) is promoted to U {Ni\N2), and they can be defined 
in any super-representation TZ. In [16] it has been argued that this 1/2 BPS loop, which we will 
denote by SW-ji, localizes to the matrix model correlator 

ABJM 



sn\ V -e"^' 

in the ABJM matrix model. Here, 



sn = Stm (^J . (2.11) 

When 7?. = □, we have the simple relationship 

(ST^n) = ^^1^^ (iVi(W^n) + N2{Wa)) . (2.12) 

In general, as it is clear from (2.10), the vev of the 1/2 BPS Wilson loop can be obtained if we 
know the vevs of the 1/6 BPS Wilson loop, but we expect it to be simpler. 

The work of [5] reduced the computation of vevs of Wilson loops in ABJM theory to the 
computation of matrix model correlators in the matrix model (2.6). Perturbative calculations 
are now straightforward. However, in order to obtain exact interpolation functions, we have to 
resum all double-line diagrams at fixed genus in the matrix model. This is straightforward in 
the Gaussian matrix model which computes the vev of the 1/2 BPS Wilson loop in A/" = 4 SYM 
[3, 2, 1], but it is not for the model (2.6). As we will see in this paper, the most efficient way to 
solve this matrix model in the l/N expansion is to relate it to the lens space matrix model of 
[12] and to its large A'^ string dual [13]. 

3. The lens space matrix model and its large N dual 

The lens space matrix model was introduced in [12, 13] in order to compute the partition function 
of Chern-Simons theory on lens spaces of the form L{p, 1) = S^/Zp. In this paper we will be 
particularly interested in the case p = 2. In this case, the matrix model has the structure 



(3.1) 



X 

'■,J 
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and it describes the expansion of the Chern-Simons partition function around a generic non- 
trivial flat connection, corresponding to the symmetry breaking pattern 



U{N) U{Ni) X UiN2). 
The model has a large expansion of the form 

oo 

F = logZ = J2F9{ti,t2)9!'-^ 

9=0 

where 

ti = QsNi 

are the partial 't Hooft parameters. The genus zero free energy has the structure 
where F^{t) is the Gaussian matrix model genus zero amplitude, 



1 



Fo^(t) = -t-(logt 



(3.2) 

(3.3) 

(3.4) 
(3.5) 

(3.6) 



and (ti,t2) is the contribution from fatgraphs of genus zero. The first nontrivial terms are 



F^{h,t2) = —(tf + 6tlt2 + 18tl4 + 6h4 + 4) 



1 



(3.7) 



{4t1 + 45*5**2 + 225tftl + 1500^4 + 225t^t^ + 45tit^ + 4^) + • 



3^.3 
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Higher genus free energies can be computed analogously. 
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Figure 1: The cuts for the CS lens space matrix model in the z plane and in the Z — plane. 

The matrix model (3.1) can be studied from the point of view of the 1 /N expansion. This was 
done in detail in [13], where this expansion was identified with the genus expansion of topological 
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string theory on a local Calabai-Yau manifold. Many results of [13] were rederived in the paper 
[14, 19] by using standard matrix model techniques, which we now review. 

At large N, the two sets of eigenvalues, Hi, uj, condense around two cuts centered around 
z = 0, z = iri, respectively. We will write them as 

Ci = {-A,A), C2 = (7ri-S,7ri + 5), (3.8) 

in terms of the endpoints A, B. It is also useful to use the exponentiated variable 

Z = e^ (3.9) 

In the Z plane the cuts (3.8) get mapped to 

(l/a,a), (-1/6,-6), (3.10) 

which are centered around Z = 1, Z = — 1, respectively, and 

a = e^, 6 = e^, (3.11) 

see Fig. 1. We will use the same notation Ci,2 for the cuts in the Z plane. An important quantity 
introduced in [14] is the total resolvent of the matrix model, oj{z). It is defined as 

^{^) = 9s ^Ecoth ^i+gs (f:tanh ^ . (3.12) 

In terms of the Z variable, it is given by 

a;(.)d. = -(ii + t2)^ + 2gs + 2gs (3.13) 

and it has the following expansion as Z — > 00 

u;{z)^h+t2 + ^(j2^''^-J2^''n + --- (3-14) 

\ i=l j=l I 

From the total resolvent it is possible to obtain the density of eigenvalues at the cuts. In the 
planar approximation, we have that 

u;o{z) = -ih+t2) + 2ti [ pi(M)^^d/x + 2t2 / P2(z^)^^di., (3.15) 

where pi(/x), p2{i^) are the densities of eigenvalues on the cuts Ci, C2, respectively, normalized as 

/ Pi(//)dM= / p2i'^)diy = l. (3.16) 

JCi JC2 
The standard discontinuity argument tells us that 

1 d.X 

pi{X)dX = --r^^ (wo(X + ie) - MX - ie)) , X G Ci, 

(3.17) 

P2(l^)dy= -A-^(a;o(^ + ie)-u;o(i^-ie)), Y eC^,. 
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The planar resolvent (3.15) was found in explicit form in [14]. It reads, 



u!o{z) = 2 log 



e-*/2 



where 



t = ti+t2 

is the total 't Hooft parameter. It is useful to introduce the variables 

a = a + -, (3 = b+Y. 
a 



as well as 



c = 



a- 13 



a + P 



2 ' ^ 2 

The expansion (3.14) implies then 

C = 2e*. 

As it is standard in matrix models, the 't Hooft parameters turn out to be period integrals with 
a nontrivial relation to a, b: 



(3.18) 
(3.19) 

(3.20) 

(3.21) 
(3.22) 



= <^ u)o{z)dz, t2 = - — r (f uJo{z)dz. 



47ri Jc 



Ci 



47ri Jc 



C2 



(3.23) 



The derivatives of these periods can be calculated in closed form by adapting a trick from [20]. 
If we write 



UJq{z) = log 



with 

it follows that 



f{z) = z^-Cz + i, 



1 



dZ 



ab 



dC 47ri 7ci,2 y^fiZy-eZ^ 7r(l + ab) 

where K{k) is the complete elliptic integral of the first kind, and its modulus is given by 

,_ (a^-l)(6^-l) 
(1 + a6)2 ■ 



(3.24) 
(3.25) 
(3.26) 

(3.27) 



The above relationships determine in principle the planar content of the theory. However, 
this matrix model solution is further clarified by considering its large dual, which was in 
fact discovered before [13]. This dual is given by topological string theory on the anti-canonical 
bundle of the Hirzebruch surface Fq = P-*^ x P-*^. The mirror geometry is encoded in a family of 
elliptic curves S, which can be written as 



y 



z\x^ + a; + 1 — ^(1 + X + z\x^Y ~ ^z^x^ 



(3.28) 



Here, zi,Z2 parametrize the moduli space of complex structures, which is the mirror to the 
enlarged Kahler moduli space of local Fq. This moduli space has a very rich structure discussed 
in [13, 21]. 
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The mirror geometry (3.28) is nothing but the spectral curve of the lens space matrix model, 
and it is closely related to the resolvent uJo{Z). Indeed, one finds that loo(Z) ~ log y(x) provided 
we identify the variables as 

X = -Zzl^l"^, (3.29) 



C = i = iJ^. (3.30) 



and 



In order to make further contact with the matrix model, one has to look at the moduli space 
of (3.28) near the orbifold point discovered in [13]. This is defined as the point xi = X2 = in 
terms of the variables: 

-21 1 
xi = 1 , X2 = 7 ^. (3.31) 

V^(l-l) 

Using mirror symmetry, we can calculate the periods of uJo{z) along the cycles of the spectral 
curve as solutions to a Picard-Fuchs equation. In terms of the coordinates xi,2, the Picard-Fuchs 
system is given by the two operators 

£1 = 1 (8 - 8x1 + xj) X2d:c, + I (-4 + (-2 + xi) \i) dl^ 

+ (-1 + xi) xfdxi - xi (2 - 3x1 + xf) X2dx-^dx2 + (-1 + xi) ^xjdl^, 

£2 = (2 - xi) X2dx2 + (-1 + (1 - xi) x|) - xjdxi 
+ 2 (-1 + xi) xiX2dxidx2 + (1 - xi) x\dl^. 

A basis of periods near the orbifold point was found in [13]. It reads, 

(Ti = - log(l - Xi) = ^ CmfiX^, 
m 

C2 = ^ ^ (hn,nX'ix'^, 33^ 

m,n 

^(T2 = log(xi) + ^ dm,nXfx2 , 

where the coefficients Cm,n S'lid dm,n are determined by the following recursions relations 

(n + 2 - 2mY 



^m,n — I'm— l,n 



4(m — n)(m — 1) ' 

(Cm,n-2(ra - m - l)(n - m - 2) - Cra-l,n-2{n - m - 1)^), 



n(n — 1) 

_d^_i,„(n + 2 - 2rrif + 4(77, + 1 - 2m)cm,n + 4(2m - n - 2) (3.34) 

4(m — n)(m — 1) 

dm,n = , ^ -,M m,n-2{n -m - l){n - Ul - 2) - dm-l,n-2in - Ul - if 

n(n — 1) 

+ (2n - 2 - 2m)cm_i,n-2 + (2m + 3 - 2n)cm,n-2)- 
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The 't Hooft parameters of the matrix model are related to the periods above as 

ii = ^(cji + (J2), t2 =^{(Ji-a2). (3.35) 

The remaining period in (3.33) might be used to compute the genus zero free energy of the matrix 
model. Notice that xi^2 (or equivalently C;^ defined in (3.21)) are "bare" coordinates, while 
(71^2 (and therefore ti^2) are flat coordinates, annihilated by the Picard-Fuchs operators. 

It is now a matter of (computer) routine to calculate the different quantities, like the end- 
points of the cuts, as an expansion in the 't Hooft parameters. One obtains, for example, 

a = 1 + 2^/h. + 2ti + ^^^1(3*1 + ^2) + ti{ti + t2) + --- (3.36) 
The expansion for b is obtained from this one simply by exchanging ti 0*2- 



4. Exact results for the ABJM model 



4.1 General results 



We will now use our knowledge of the solution of the lens space matrix model to solve the ABJM 
model, at least at the planar level. It is clear that the matrix model (2.6) is closely related to 
(3.1), but there are some obvious differences: in (2.6) the interaction between the /x and the 
u eigenvalues is in the denominator, and the Gaussian action for the vs has the opposite sign. 
These ingredients are precisely the ones needed to make (2.6) a supergroup extension of (3.1). We 
will now quickly review some results on supermatrix models, following [22, 23, 24]. A Hermitian 
supermatrix has the form 



where A (C) are iVi x A^^i (A^2 x -^2) Hermitian, Grassmann even matrices, and ^ is a complex, 
Grassmann odd matrix. The supermatrix model is defined by the partition function 

Z^{Ni\N2) = J P$e"^^*'^^*^ (4.2) 

where we consider a polynomial potential V{^), and Str is the supertrace 

Str $ = Tr A - Tr a (4.3) 

There are two types of supermatrix models with supergroup symmetry U(Ni\N2). the ordinary 
supermatrix model, and the physical supermatrix model [23]. The ordinary supermatrix model is 
obtained by requiring A, C to be real Hermitian matrices, while the physical model is obtained by 
requiring that, after diagonalizing $ by a superunitary transformation, the resulting eigenvalues 
are real. Here we will be interested in the physical supermatrix model. Its partition function 
reads, in terms of eigenvalues [23, 24] 

z.{N,\N2) = f n n d^.- " ^" - ^^''^^ • (4-4) 
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When the two groups of eigenvalues yUj, Vj are expanded around two different critical points, the 
partition function (4.4) is well-defined as an asymptotic expansion in Qs- It is easy to show that 
(4.4) is related to the partition function of the corresponding bosonic, two-cut matrix model 

z^{N,,N,) = / n d/.. n d^.- n - ^'^f - n (/^^ - "^^''^^^^ ""^^'^^ 

i=\ j=l i<j i,j 

(4.5) 

after changing N2 —N2: 

Z,{Ni\N2) = Z^{Ni,-N2). (4.6) 

Such a flip of sign is trivially performed if one knows the exact solution of the model in the 
1/iV expansion. The relation (4.6) can be proved diagramatically by introducing Faddeev-Popov 
ghosts as in [24, 25]. 

We now see that the relationship between the ABJM matrix model and the lens space matrix 
model is identical to the one we have between supergroup matrix models and multi-cut bosonic 
matrix models, with the only difference that the interaction between the eigenvalues has been 
promoted to the sinh interaction typical of Chern-Simons matrix models. Indeed, the lens space 
matrix model is a two-cut matrix model where the /x, u eigenvalues are expanded around two 
different saddle points, z = and z = it\. The ABJM matrix model is just its supergroup version. 
We then conclude that 

ZABm{Ni,N2,g) = Zn2,i){Ni,-N2,g). (4.7) 
Notice that the change N2 ^ —N2 is equivalent to setting 

ti = 2mXi, t2 = -2mX2, (4.8) 

where 

A, = i = l,2, (4.9) 

are the 't Hooft parameters of the ABJM model. 

The appearance of a hidden supergroup structure in the matrix model of [5] is not surprising, 
since A/" = 4 Chern-Simons-matter theories are classified by supergroups [26]. In fact, the ABJM 
theory can be constructed as an A^ = 4 theory with supergroup U{Ni\N2) and containing both 
hypermultiplcts and twisted hypermultiplets [27] . This hidden supergroup structure in the ABJM 
theory is explicitly used in the construction of half-BPS Wilson loops in [16]. 

Let us now discuss Wilson loops. The most natural correlator in the standard lens space 
matrix model is 



-0"^ , 

where 7?. is a representation of U{Ni + N2). In Chern-Simons gauge theory on L(2, 1), this 
computes the vev of a trivial knot, expanded around a generic, fixed flat connection. In the 
topological string large N dual, it corresponds to an open string amplitude for a toric D-brane 
(see for example [18] for more details). These vevs can be computed for any TZ and to all orders 
in the 1/N expansion [17, 18]. 

In order to consider its supergroup extension, notice that a representation of U{Ni + A^2) 
induces a super-representation of U{Ni\N2), defined by the same Young tableau TZ (see for 
example [28]). Therefore, the supergroup generalization of (4.10) is simply 

(^'- Co 
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This can be also written as [28] 



k i \ 



= E — n(^ {^''') - (-1)'^ («'"^^'' 



(4.12) 



In this equation, which is the supergroup generahzation of Probenius formula, A; = (ki) is a vector 
of non-negative, integer entries, which can be regarded as a conjugacy class of the symmetric 
group, XTl{^) is the character of this conjugacy class in the representation TZ, and 



^k 



Yli'''ki\ (4.13) 

e 

We then have 

which extends (4.7) to correlation functions. In view of (2.10), we conclude that the vevs of the 
1/2 BPS Wilson loops constructed in [16] can be computed to all orders in the 1/N expansion 
by calculating the correlator (4.10) in the lens space matrix model and changing N2 — > —N2 

{SWtz) = - (Trn _l,^)){N^,-N2,g). (4.15) 



SR 



Let us give an example of how to calculate the vev (4.15) when TZ = n. In the planar limit 
we have, for the lens space matrix model correlator (4.10), the exact answer 

' /T^off ("6, 



A^i - iV2 \ V -e'^V / ii - t2 2 ' 
where ( is defined in (3.21). After setting (4.8) we obtain the weak coupling expansion 

{SWu) = l+i7r(Ai-A2)- Jtt^ (2X1 " ^AsAi + 2Xl)-\m^ (A? - 4A2A? + 4A^Ai - Ai)+- • ■ (4.17) 

This is of course (up to normahzations) the first term in equation (6.54) of [18]. Higher genus 
corrections can be extracted from the higher genus resolvents ojg{z). These in turn can be 
computed by using standard matrix model techniques applied to the spectral curve (3.28), as in 
[17, 18]. 

4.2 The 1/6 BPS Wilson loop 

If the densities of eigenvalues pi{p) and p2(i^) given in (3.17) are known, it is possible to calculate 
the exact planar limit of the correlator 

^(IVe'^^)= f pi(Ai)e'*d/x= f piiX)XdX, (4.18) 
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as well as of multiple-winding correlators 

l(IVe"'^0= / Pi(/x)e"^dM. 



(4.19) 



We then conclude, in view of (2.8), that the planar limit of the vev of a 1/6 BPS Wilson loop is 
given by 

(^□) = / Pi(M)e'^dM (4.20) 

after changing variables as in (4.8). 

The densities pi(/i) and p^iy) can be explicitly calculated from (3.17) and (3.18). We find. 



px[X)^X = tan-i 

irti 



P2{Y)dY = tan-i 

In terms of the variable x = log X we have 

1 

pi{x) = tan 



px + i + x^ 



(3Y + 1 + Y^ 

ay - 1 - y2 



dX 

dY 
~Y' 



(4.21) 



/ a — 2 cosh X 
(3 + 2 cosh X 



(4.22) 



and a similar expression for P2{y)- Notice that, if t2 = 0, one has /3 = 2, a = 4e* — 2, and pi{x) 
becomes the density of eigenvalues for the matrix model of Chern-Simons theory on [29] 



e* - cosh^ f 5' 



P-^-) = Vt'''''' cosh (I) 

Remark 4.1. We can write the density of eigenvalues (4.22) as 

Pi{x) = f{x)-\/A^ - x2, 

where 

fix) = , tan" 

We find for example 



irtiV A'^ — X' 



I a — 2 cosh X 
/3 + 2 cosh X 



k=0 



X 



2k 



ttUA 



tan 



a -2 
J+2 



(4.23) 

(4.24) 
(4.25) 

(4.26) 



It is easy to check that pi{x) agrees, in the special case ti = —t2 = t, with the perturbative expan- 
sion obtained in [11] up to order 10 in t. Our expression (4.22) gives then the full resummation 
of the expression obtained in [11], and extends it to any ti,t2- 



The integral (4.20) is then given by 
{Wu) = ^h, h-- 



na 

I tan"^ 

Jlla 



laX -1-^2 

px + i + x^ 



dX. 



(4.27) 
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This integral is not easy to calculate in closed form, but its derivatives w.r.t. and ^ can be 
expressed in terms of elliptic integrals. We find 



^ ^ ^"^^ ^ {aK{k)-{a + b)U{n\k)). 



dC 27 ^(aX - 1 - X^){I3X + 1 + X2) V^{l + ab) 

l/a 

dh_lf + l + X^)-{aX -1- X^)) dX _ Va6 

l/a 

(4.28) 

where n(n|fc) is the complete elliptic integral of the third kind, K{k),E{k) are elliptic integrals 
of the first and second kind, respectively, the modulus is given by (3.27), and finally 

ha? -I 

n = -- -. 4.29 

al + ah ^ ' 

This determines the planar limit of the 1/6 BPS Wilson loop exactly. As we mentioned in the 
introduction, the Wilson loop is naturally expressed in terms of the "bare" coordinates and 
we have to use the mirror map (3.33) in order to write it in terms of the 't Hooft parameters. 

As an application of these formulae, wc will present the first few terms of the weak coupling 
expansion in ti,t2- To do this, we simply calculate 

du ~ dCdti d^^' ^ '^^^ 

we use the expressions (4.28), and we integrate w.r.t. ti, t2- This produces the expansion 

( Wq) = 1 + 1 + ^ {2tl + 3*1*2) + ^ (2t? + 6*2*2 + 4*1*2) + ^ (8*f + 35*?*2 + 30*?*^ + 10*i*i) + " " " 

(4.31) 

The result for the Wilson loop in the ABJ theory is obtained by simply changing variables to 
(4.8). The result agrees at the first few orders with a perturbative matrix model calculation of 
[16]. Finally, the vev of the other Wilson loop (2.9) is obtained by exchanging *i 0*2, or, in the 
ABJ theory, by exchanging Ai ■<->■ A2 and complex conjugating the result. 

Interestingly, the strong coupling limit of the above expressions depends on the direction in 
which one goes to infinity. We will discuss one such direction in detail in the next subsection, 
where we consider the restriction to the original ABJM model A^i = N2. 

4.3 The ABJM slice 

In the original ABJM theory with Ni = N2 = N (the case A''i / A^2 was considered in [7]) we 
should look at the slice 

N 

*i = -*2 = 27riA, A = — (4.32) 

k 

in the moduli space of the dual topological string. From the point of view of the periods ai , a"2 
in (3.33) this means that we should set 

(Ti = 0, (4.33) 
therefore xi = 0. In order to have a nontrivial (72, we must consider the double-scaling limit 

xi 0, a;iX2 = C fixed. (4.34) 
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The one- dimensional subspace (4.32) corresponds, in terms of the variables ^, to ^ = 2. As in 
[30], we can find simplified expressions for the periods in this subspace. It is easy to see from the 
structure of (72 that, in the limit (4.34), one has 

oo 

0'2 = ^2 ^rnC^^'^^, am = C2m+l,2m+l, (4.35) 
m=0 

and from the recursion relation (3.34) we find 

am = ^ 4.36 

7r(2m + l)r (m + 1)2 ^ ^ 

We then obtain 

which is in fact a particular case of (3.26), as it can be easily seen by using the transformation 
properties of the elliptic integral K{k). The period ti itself can be written as a generalized 
hypergeometric function: 

^i(C) = 7 3i^2fiU;i>?;SV (4-38) 



4 " V2 2' 2' ' 2' 16, 

In the physical ABJM theory, ti is purely imaginary. This means that ^ is purely imaginary as 
well, so we set 

C = iK (4.39) 



and we finally obtain 



Let us now calculate the planar limit of the vev of the 1/6 BPS Wilson loop in the ABJM 
slice Ni = N2. Since 

a = 2 + iK, ^ = 2 - iK, (4.41) 
the endpoints of the cuts are given by 

a{K) = ^ (2 + Ik + 1/^(41- k)) , 

^ ' (4.42) 

6(k) = - (^2 - Ik + i/-«(4i + «)) • 

The planar vev of the Wilson loop is then determined, as a function of the 't Hooft coupling A, 
by the single equation 

dK 277^ Vab (l + ab) 

together with the explicit relation between A and k in (4.40) -yet another example of mirror 
map. 

As a check, we can perform a weak coupling expansion. The weakly coupled region corre- 
sponds to 

K < 1, A < 1, (4.44) 
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Figure 2: Plot of (4.52) as a function of A, displaying clearly an interpolation between a quadratic 
behavior near the origin and a square-root growth as A becomes large. 

and in this region the variables are related as 

8^-^+^ 6^ + ^260 18144^ + ^^^ ^' ^^'^^^ 

which is obtained from the inversion of (4.40). By expanding (4.43) in power series around k = 0, 
and using the mirror map (4.45), we obtain 

^ °' V 6 2 120 12 1008 10 362880 

2603i7r9A9 344739l7riOAio 1166161 ivriUii 52393723197ri2Ai2 ,^,,...\ 
15120 7983360 3628800 6227020800 ^ V 

(4.46) 

where we have extracted a framing factor to facilitate comparison with existing results like (2.4). 
The first few terms agree with the calculations in [8, 9, 5, 11]. 

Of course, the main advantage of having analytic expressions is that one can perform a 
weak-strong coupling interpolation easily. The strong coupling region is 

K > 1, A > 1 (4.47) 

and (4.40) leads to the asymptotic expansion 

which is immediately inverted to 

^ = e-v^ (l + o( ^,e-2-v^^^ . (4.49) 



.Va 

On the other hand, it is easy to check from (4.43) that 
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{Wu)--^e-^\ A»l. (4.51) 



It follows that 



The leading exponential is in perfect agreement with the AdS prediction (2.5), and the exact 
answer interpolates between the weak and the strong coupling behaviors. This is illustrated in 
Fig. 2, which represents the function 

f{\) = \og\{Wn)\. (4.52) 

This function interpolates between 

^ 57r2A2 



-, A<1, 

(4.53) 



/(A) = < 

^7r\/2A, A»l. 

Let us now consider the vev of the 1/2 BPS Wilson loop of [16] in the fundamental repre- 
sentation, which is given by (2.12). The expression for this vev follows from the specialization of 
(4.16), 

Since, in the ABJM slice, (2.8) and (2.9) are related by complex conjugation, this vev is manifestly 
real. Notice also that, when adding the contributions of the two gauge groups as in (2.12), all 

terms in the series expansion at small k of the r.h.s. of (4.43) cancel except for the first constant 
term. This confirms that the 1/2 BPS Wilson loop is much simpler than the 1/6 BPS Wilson 
loop. At strong coupling we have, 

(St^n)~^e-^, A»l, (4.55) 

which displays the same leading exponential behavior predicted by the large A*" dual. 

The computations at strong coupling can be easily extended to the case Ni ^ N2, but they 
depend on the direction in which we take the limit in the space of 't Hooft parameters. For ^ 
fixed and C, large (therefore Ai A2), we find the same exponential behavior 



5. Conclusions 



In this paper we have related the matrix integral of [5], computing Wilson loop vevs in ABJM 
theory, to a supergroup extension of the lens space matrix model introduced in [12] and solved in 
[13, 14]. This has made possible to obtain exact expressions for the planar vev of 1/6 BPS Wilson 
loops, and we have verified the strong coupling behavior predicted by the AdS dual [8, 10, 9]. 

There are various avenues for further research. It is clearly gratifying to see that the tech- 
niques of mirror symmetry are relevant to the strong- weak interpolation problem in the AdS/CFT 
correspondence, but it would be nice to have an a priori understanding of the relationship be- 
tween topological string theory and the ABJM theories, maybe along the lines of [31]. This would 
lead to further fruitful interactions between the two topics. 
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We should also mention that the strong coupling analysis of Wilson loops in the context 
of the gauge theory is different from what one does in topological string theory. There, the 
interpolation from weak to strong coupling involves a different choice of flat coordinates and 
a different choice of duality frame (see for example [32] and specially the analysis of Wilson 
loop amplitudes in [18]). This means, in particular, that in topological string theory, analytic 
continuation of the amplitudes is not enough. In the gauge theory analysis, in contrast, one 
uses the same flat coordinates (the 't Hooft parameters) and the same duality frame in the full 
moduli space, and the weak and the strong coupling regions are simply related by an analytic 
continuation. 

One obvious problem is to find expressions for the Wilson loop vevs beyond the planar 
approximation. As we have explained in this paper, this can be trivially done for the 1/2 BPS 
Wilson loops constructed in [16], by simply applying the techniques of [17, 18]. For the 1/6 BPS 
Wilson loop the calculation of 1/A^ corrections is more difficult, but in principle it can be done. 
Finally, it would be interesting to compare the matrix model results with an AdS calculation 
using D-branes, as it was done in [33] for the 1/2 BPS Wilson loop of A/" = 4 Yang-Mills theory. 

Acknowledgments 

We would like to thank Anton Kapustin for discussions. We are specially grateful to Nadav 
Drukker for many conversations which helped us very much, and we want to thank him, as well 
as Diego Trancanelli, for sharing their unpublished results with us. This work is supported in 
part by the Fonds National Suisse. 

References 

[1] V. Pestun, "Localization of gauge theory on a four-sphere and supersymmetric Wilson loops," 
arXiv:0712.2824 [hep-th]. 

[2] N. Drukker and D. J. Gross, "An exact prediction of N = 4 SUSYM theory for string theory," J. 
Math. Phys. 42, 2896 (2001) [arXiv:hep-th/0010274]. 

[3] J. K. Erickson, G. W. Semenoff and K. Zarembo, "Wilson loops in N = 4 supersymmetric 
Yang-Mills theory," Nucl. Phys. B 582, 155 (2000) [arXiv:hep-th/0003055]. 

[4] G. W. Semenoff and K. Zarembo, "Wilson loops in SYM theory: From weak to strong coupling," 
Nucl. Phys. Proc. Suppl. 108, 106 (2002) [arXiv:hcp-th/0202156]. 

[5] A. Kapustin, B. Willett and 1. Yaakov, "Exact Results for Wilson Loops in Superconformal 
Chern-Simons Theories with Matter," arXiv:0909.4559 [hep-th]. 

[6] O. Aharony, O. Bergman, D. L. Jafferis and J. Maldacena, "N=6 superconformal 

Chcrn-Simons-mattcr theories, M2-branes and their gravity duals," JHEP 0810, 091 (2008) 
[arXiv:0806.1218 [hep-th]]. 

[7] O. Aharony, O. Bergman and D. L. Jafferis, "Fractional M2-branes," JHEP 0811, 043 (2008) 
[arXiv:0807.4924 [hep-th]]. 

[8] N. Drukker, J. Pleflsa and D. Young, "Wilson loops in 3-dimensional N=6 supersymmetric 

Chern-Simons Theory and their string theory duals," JHEP 0811, 019 (2008) [arXiv:0809.2787 
[hep-th]]. 

[9] B. Chen and J. B. Wu, "Supersymmetric Wilson Loops in N=6 Super Chern-Simons-matter 
theory," Nucl. Phys. B 825, 38 (2010) [arXiv:0809.2863 [hep-th]]. 



-17- 



[10] S. J. Rey, T. Suyama and S. Yamaguchi, "Wilson Loops in Superconformal Chcrn-Simons Theory 
and Fundamental Strings in Anti-de Sitter Supergravity Dual," JHEP 0903, 127 (2009) 
[arXiv:0809.3786 [hep-th]]. 

11] T. Suyama, "On Large N Solution of ABJM Theory," arXiv:0912.1084 [hep-th]. 

12] M. Marifio, "Chcrn-Simons theory, matrix integrals, and perturbative three-manifold invariants," 
Commun. Math. Phys. 253, 25 (2004) [arXiv:hep-th/0207096]. 

13] M. Aganagic, A. Klemni, M. Marino and C. Vafa, "Matrix model as a mirror of Chern-Simons 
theory," JHEP 0402, 010 (2004) [arXiv:hep-th/0211098]. 

14] N. Halmagyi and V. Yasnov, "The spectral curve of the lens space matrix model," JHEP 0911, 104 

(2009) [arXiv:hep-th/0311117]. 

15] R. Gopakumar and C. Vafa, "On the gauge theory /geometry correspondence," Adv. Theor. Math. 
Phys. 3, 1415 (1999) [arXiv:hep-th/9811131]. 

16] N. Drukker and D. Trancanelli, "A supermatrix model for = 6 super Chern-Simons-matter 
theory," arXiv:0912.3006 [hep-th]. 

17] M. Marino, "Open string amplitudes and large order behavior in topological string theory," JHEP 
0803, 060 (2008) [arXiv:hep-th/0612127]. 

18] V. Bouchard, A. Klemm, M. Marino and S. Pasquetti, "Remodeling the B-model," Commun. Math. 
Phys. 287, 117 (2009) [arXiv:0709.1453 [hep-th]]. 

19] N. Halmagyi, T. Okuda and V. Yasnov. "Large N duality, lens spaces and the Chern-Simons matrix 
model," JHEP 0404, 014 (2004) [arXiv:hep-th/0312145]. 

20] A. Brini and A. Tanzini, "Exact results for topological strings on resolved Y(p,q) singularities," 
Commun. Math. Phys. 289, 205 (2009) [arXiv:0804.2598 [hep-th]]. 

21] B. Haghighat, A. Klemm and M. Rauch, "Intcgrability of the holomorphic anomaly equations," 
JHEP 0810, 097 (2008) [arXiv:0809.1674 [hep-th]]. 

22] L. Alvarez-Gaume and J. L. Maiies, "Supermatrix models," Mod. Phys. Lett. A 6, 2039 (1991). 

23] S. A. Yost, "Supermatrix models," Int. J. Mod. Phys. A 7, 6105 (1992) [arXiv:hep-th/91 11033]. 

24] R. Dijkgraaf and C. Vafa, "N = 1 supersymmetry, deconstruct ion, and bosonic gauge theories," 
arXiv:hep-th/0302011. 

25] R. Dijkgraaf, S. Gukov. V. A. Kazakov and C. Vafa, "Perturbative analysis of gauged matrix 
models," Phys. Rev. D 68, 045007 (2003) [arXiv;hep-th/0210238]. 

26] D. Gaiotto and E. Witten, "Janus Configurations, Chern-Simons Couplings, And The Theta-Angle 
in N=4 Super Yang-Mills Theory," arXiv:0804.2907 [hep-th]. 

27] K. Hosomichi, K. M. Lee, S. Lee, S. Lee and J. Park, "N=5,6 Superconformal Chern-Simons 
Theories and M2-branes on Orbifolds," JHEP 0809, 002 (2008) [arXiv:0806.4977 [hep-th]]. 

28] I. Bars, "Supergroups And Their Representations," Lectures Appl. Math. 21, 17 (1983). 

29] M. Marino, "Les Houches lectures on matrix models and topological strings," arXiv:hep-th/0410165. 

30] M. Marifio, S. Pasquetti and P. Putrov, "Large N duality beyond the genus expansion," 
arXiv:0911.4692 [hep-th]. 

31] G. Bonclli and H. Safaai, "On gauge/string correspondence and mirror symmetry," JHEP 0806, 
050 (2008) [arXiv:0804.2629 [hep-th]]. 

32] M. Aganagic, V. Bouchard and A. Klemm, "Topological Strings and (Almost) Modular Forms," 
Commun. Math. Phys. 277, 771 (2008) [arXiv:hep-th/0607100]. 

33] N. Drukker and B. Fiol, "All-genus calculation of Wilson loops using D-branes," JHEP 0502, 010 
(2005) [arXiv:hep-th/0501109]. 



- 18 - 



